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We present a linearized shear-transformation-zone (STZ) theory of glassy dynamics in which the 
internal STZ transition rates are characterized by a broad distribution of activation barriers. For 
slowly aging or fully aged systems, the main features of the barrier-height distribution are deter- 
mined by the effective temperature and other near-equilibrium properties of the configurational 
degrees of freedom. Our theory accounts for the wide range of relaxation rates observed in both 
structural glasses and soft glassy materials such as colloidal suspensions. We find that the frequency 
dependent loss modulus is not just a superposition of Maxwell modes. Rather, it exhibits an a peak 
that rises near the viscous relaxation rate and, for nearly jammed, glassy systems, extends to much 
higher frequencies in accord with experimental observations. We also use this theory to compute 
strain recovery following a period of large, persistent deformation and then abrupt unloading. We 
find that strain recovery is determined in part by the initial barrier-height distribution, but that true 
structural aging also occurs during this process and determines the system's response to subsequent 
perturbations. In particular, we find by comparison with experimental data that the initial defor- 
mation produces a highly disordered state with a large population of low activation barriers, and 
that this state relaxes quickly toward one in which the distribution is dominated by the high barri- 
ers predicted by the near-equilibrium analysis. The nonequilibrium dynamics of the barrier-height 
distribution is the most important of the issues raised and left unresolved in this paper. 
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I. INTRODUCTION 

A growing body of experimental evidence indicates 
that glassy materials exhibit broad spectra of linear re- 
laxation modes. Most definitively, the frequency depen- 
dent loss modulus, G"{uj), exhibits a broad peak that 
has a maximum near the viscous relaxation rate and ex- 
tends over two or more decades of higher frequencies. 
This behavior has been seen in structural and metallic 
glasses 0-9, and also in soft glassy materials such as 
colloidal suspensions P-fllj. Additional evidence for a 
broad relaxation spectrum emerges from strain recovery 
experiments, in which materials are unloaded abruptly 
after shear deformation d, [l^, [111 ■ 

Our intention in this paper is to develop a theory that 
quantitatively accounts for these measurements, extend- 
ing the results announced in [3]. To date, two of the 
most successful descriptions of these phenomena have 
been the soft glassy rheology (SGR) theory of Sollich, 
Gates, et al. [15|, [la], and the attempt by Gates and 
coworkers to use mode-coupling theory (MGT) for simi- 
lar purposes [13, [Si- Here, we base our analysis on the 
shear-transformation-zone (STZ) theory of amorphous 
plasticity 19-.22i]. We contend that the STZ theory pro- 
vides a more suitable starting point than either SGR or 
MGT. Ultimately, however, our goal is to find a descrip- 
tion of linear glassy dynamics that combines the strongest 
elements of all three of these theories. 

In contrast to SGR, which postulates a trap-like mech- 
anism, the STZ theory is based directly on a model of the 
localized shear transformations that are observed in col- 
loidal suspensions, bubble rafts, and molecular dynam- 



ics simulations, and that presumably occur in molecu- 
lar glasses as well. The STZ theory also describes the 
nonequilibrium flow of energy and entropy in a way that 
has not yet been achieved in SGR. On the other hand, 
the distribution of trapping energies that is an essential 
component of SGR has, until now, been missing in the 
STZ theory. Incorporating this feature of SGR into the 
STZ theory is a main theme of the present paper. 

In contrast to MGT, which starts with a fluidlike 
many-body Hamiltonian and makes decoupling approx- 
imations, the STZ theory is based on a solidlike pic- 
ture of an amorphous system. This is an intrinsically 
non-perturbative strategy. Among other phenomena, 
this strategy allows the STZ theory to predict a sharp 
dynamic transition between jammed and flowing states 

- a transition that seems difficult to achieve in any 
perturbation-theoretic analysis that starts with liquid- 
like variables - but which is accomplished at least in part 
in [l3, [3 . The success of the STZ theory comes at the 
expense of a seemingly irreducible degree of phenomenol- 
ogy. Without a precise basis in many-body physics, the 
STZ theory (like SGR) must postulate a set of relevant 
internal dynamic variables and use whatever constraints 
are available - symmetry principles, the laws of thermo- 
dynamics, physical insight, and agreement with experi- 
ment - to deduce equations of motion for those variables. 

The STZ theory is based on two fundamental assump- 
tions. The first of these is that the degrees of freedom of 
a glassy material can be separated into two subsystems 

- the slow configurational degrees of freedom, i.e. the in- 
herent structures [H, [24| , and the fast kinetic- vibrational 
degrees of freedom which, in the case of colloids, in- 
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elude the motions of the fluid in which the particles are 
suspended |25!i28j]. When external perturbations drive 
the material away from equilibrium, the effective tem- 
perature of the configurational subsystem may depart 
from the temperature of the kinetic-vibrational subsys- 
tem. The latter temperature generally is the same as the 
temperature of the thermal reservoir. 

The second fundamental premise of the theory is that 
irreversible shear deformations occur only at rare, local- 
ized, two-state, flow defects. These are the STZ's which, 
by definition, belong to the configurational subsystem. 
In flowing states, the STZ's appear and disappear as 
the driven system makes transitions between its inherent 
structures. In jammed states, the STZ's are configura- 
tionally frozen, but they still are able to make transitions 
between their internal states in response to ordinary ther- 
mal fluctuations and applied stresses. The two-state na- 
ture of the STZ's is chiefly responsible for the dynamic 
transition between jammed and flowing states at large 
driving stresses. As will be seen here in the linear regime, 
this two-state dynamics - plus a fundamentally different 
kinematic starting assumption - produces a formula for 
the frequency-dependent modulus G{uj) that is not the 
same as that which emerges from SGR. 

Earlier versions of STZ theory have been based on 
the assumption that only one characteristic kind of STZ 
is needed in order to understand stress-strain relations, 
shear banding instabilities, and the like. Thus, our STZ's 
have had only one statistically significant formation en- 
ergy, and only a single energy barrier to be crossed during 
transitions between internal orientational states. Here we 
follow fl^ and change the latter assumption in order to 
understand the dynamic measurements; that is, we need 
to assume that the STZ's occur with a distribution of 
different internal barrier heights, in much the same way 
that was deduced experimentally by Argon and Kuo in 

1980 dij. 

In the sections that follow, we propose a simple form 
for an extended, multi-barrier STZ theory, following our 
earlier report We restrict our attention to the limit 
of small external stresses, because that is the limit in 
which the response of the system is most strongly lim- 
ited by activation barriers. Accordingly, this small-stress, 
linear-response limit is where we find the clearest experi- 
mental evidence for the multi-barrier picture. We outline 
the extended STZ theory in Sec. HT] In Sec. |llT] we pro- 
vide a more detailed analysis and additional support for 
the results announced in where the frequency de- 
pendent modulus G{lli) was calculated and compared to 
the data of 0, |1] , who performed oscillatory experiments 
on a wide variety of both hard structural glasses and soft 
colloidal suspensions. 

Finally, in Sec. IIVI we further extend the analysis of 
[13] by looking at the evidence for broadly distributed 
barrier heights that emerges from strain-recovery mea- 
surements and subsequent probes of partially recovered, 
i.e. partially aged, specimens. We note that the paper 
by Belyavsky et al. |13| , dealing with strain recovery in 



metallic glasses, contains a clear description of two-state 
shear transformation zones that predates our own_Jl9|. 
We focus here on experiments by Purnomo et al Q on 
colloidal suspensions, in which the strain was measured 
following a period of large, persistent deformation and 
then abrupt unloading. We find that strain recovery is 
determined in part by the initial barrier-height distribu- 
tion, but that true structural aging also occurs during 
this process and strongly determines the subsequent lin- 
ear responses. In particular, we find that the initial de- 
formation produces a highly disordered state with a large 
population of low activation barriers, and that this state 
equilibrates quickly toward one in which the distribution 
is dominated by high barriers. 

A recurring theme thoughout Sees. Ull- IIVl is the ques- 
tion of what determines the barrier-height distribution 
and how this distribution changes during aging. To dis- 
cuss this issue quantitatively, we start by considering 
fully aged or slowly aging systems in which these distribu- 
tions approach states of quasi-equilibrium determined by 
the current values of the effective temperature and other 
relevant parameters. We argue that these distributions 
are limited by the condition that STZ transition rates 
cannot be slower than the rates of spontaneous, thermally 
activated, configurational rearrangements in glassy ma- 
terials. We then make a phenomenological guess about 
the form of the barrier-height distribution for younger 
systems at higher effective temperatures, such as those 
that occur in the strain-recovery experiments. In this 
way, we are able to make some progress in understand- 
ing the physics of these phenomena; but we have not yet 
developed an equation of motion for the STZ barrier- 
height distribution analogous to the one that serves as a 
starting assumption in SGR. 

Our main conclusions are summarized in SecfVl 

II. EXTENDED STZ THEORY 
A. STZ Equations of Motion 

As in almost all earlier presentations, it is easiest and 
physically most transparent to assume that the STZ's 
are oriented only in the ± directions relative to the shear 
stress s. We lose no generality by doing this; the tensorial 
generalizations of the equations are obvious at the end of 
the analysis, but are not needed for present purposes. 

In order to describe colloidal suspensions as well as 
molecular glasses, however, we must recognize that the 
total stress s acting on the system is the sum of partial 
stresses associated with the configurational and kinetic- 
vibrational subsystems. (See |26|].) For the suspensions, 
the kinetic-vibrational partial stress is the viscous stress 
generated by hydrodynamic forces. Therefore, 

s = sc + riK -J- (2.1) 

where sc is the configurational partial stress, and 7 is the 
total shear rate, common to both the configurational and 
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kinetic-vibrational subsystems. The notation r]K* means 
that the viscosity rji^ is a time-retarded integral opera- 
tor that becomes a function of frequency after Fourier 
transformation. 



The new feature that was introduced in is that the 
STZ's are characterized by internal barrier heights, say, 
A. Let the number of ± STZ's with given A be N±{A), 
and let the total number of (coarse-grained) molecular 
sites be N. In the limit of small applied stresses, the 
master equation for N±{A) is 

TO N±{A) = R{±sc, A) N^{A) - R{tsc, A) N±{A) 



p{0) 



Neq{A) 



N±{A) 



(2.2) 



Here, R{±sc,A))/to is the rate per STZ for transitions 
between ± orientations; it is the origin of the A depen- 
dence in the theory. Only the partial stress sc appears as 
an argument of R{±sc, A)) because the STZ's are con- 
figurational defects, tq is a fundamental time scale, for 
example, a vibration period for molecular glasses or a 
Brownian diffusion time for colloidal suspensions. 

The last two terms on the right-hand side of Eq. (12. 2p 
are the rates at which STZ's are created and annihilated 
by spontaneous thermal fluctuations. We are making a 
detailed-balance approximation in which Neq{A)/2 is the 
value approached by N±{A) in steady-state equilibrium. 
p{9) is the thermal noise strength, and 9 = ksT \s the 
bath temperature in units of energy. Mechanically gen- 
erated noise, which is a prominent ingredient of earlier 
analyses [U [l^ , is second order in the applied stress and 
therefore is negligible in this linear theory. 

As usual 1231, define the internal state variables: 



A(A) = 



.(A) 



N+{A)+N.{A) 
N 

N+{A) -iV_(A) 
N+{A)+N-{Ay 



(2.3) 



According to Eq. 
variables are 



the equations of motion for these 



Tok(A)^p{e) [A,,(A)-A(A) 
where Ae,(A) = iV,,(A)/Ar; 

To m(A) = 2 C(sc, A) [r(sc. A) - m(A) 

- p{e) m(A) - ^^^^^ m(A); 



(2.4) 



A(A) 



(2.5) 



and 



C(sc,A) 

r(sc,A) 



1 r 
2 



R{sc.A)+R{-sc,A) 

R{sc,A)-R{-sc,A) 
i?(sc,A) + i?(-sc,A)- 



(2.6) 



The total rate of plastic deformation is a superposition 
of terms of the form 

Toi^P'(A) = ^ [r{sc,A)N-{A)-R(-sc,A)N+{A) 



eoA(A)C(sc,A) r(sc,A)-m(A) 



(2.7) 



where V is the volume of the system, and vq is a molec- 
ular volume that sets the size of the plastic strain in- 
crement induced by an STZ transition. We expect eq = 
N vq/V to be a number of the order of unity. 



B. Glassy Dynamics 

Most - but not all - of the linear responses of interest 
here are associated with STZ transitions between their 
internal states. The exceptions are the last two terms 
on the right-hand side of Eq. ()2.5p . which derive from the 
annihilation and creation terms in Eq. (|2.2[) . and there- 
fore do describe configurational changes. As discussed in 
earlier papers, Agg = exp (— ez /x), where x is the effec- 
tive disorder temperature in energy units, and ez is an 
STZ formation energy that, at this stage of the discus- 
sion, may be a A-dependent quantity. Thus, in Eq. (|2.4p . 
the rate at which STZ's are spontaneously created (and 
annihilated) by thermal fluctuations is proportional to 
the factor 



Acrcation(A) OC 6 ^ ^^"^^ 

To 



(2.8) 



which has the form of an activation rate, with x pl^'Y" 
ing the role of the temperature, and p(9)/tq being the 
attempt frequency. 

Equation (|2.8p embodies some of the deepest issues in 
glass physics. It expresses our contention that the glass 
transition is intrinsically a dynamic phenomenon, rather 
than a thermodynamic phase transition. Thus, the for- 
mation energy ez has the magnitude of an ordinary in- 
teraction energy; and the Boltzmann factor exp (— ez/x) 
determines the probability of finding a fluctuation of en- 
ergy ez in the configurational subsystem once it has come 
to equilibrium at temperature %. On the other hand, the 
dynamic prefactor p{9) is a super-Arrhenius function of 
the temperature, associated with the fact that the con- 
figurational rearrangements needed to form an STZ-like 
defect involve many-body fluctuations that become in- 
creasingly com plex and unlikely as the temperature de- 
creases. (See [301 for one hypothetical picture of this 
mechanism.) Accordingly, p(9) is approximately equal 
to unity at temperatures well above the nominal glass 
temperature 9g, but it becomes very small at lower tem- 
peratures. The system undergoes a dynamically sharp 
glass transition at some temperature, say, 9o < 9g, if - 
as in the Vogel-Fulcher-Tamann formula - p{9) actually 
vanishes at and below ^o- 

The difference between x ^-iid ^ is a measure of the 
extent to which the configurational degrees of freedom 
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have fallen out of equilibrium with the heat bath. In the 
absence of mechanically generated noise, configurational, 
i.e. "structural," relaxation is governed by an equation 
of motion for x, which we write in the form 



ez \ 6 



(2.9) 



Here k is a dimensionless constant, very roughly of the 
order of unity, and exp { — ca/x) is a measure of the pop- 
ulation of defects, with formation energy ca, that enable 
energy transfer from the configurational to the kinetic- 
vibrational degrees of freedom. The energy ea should be 
roughly the same as ez- Thus, Eq. (|2.9p describes slow, 
structural aging during which x relaxes toward 9. 

With these assumptions, we introduce a normalized 
distribution over barrier heights, p(A), and use Ea. (|2.7p 
to write the total rate of plastic deformation in the form 



(A) J dAp{A)C{sc,A) [T{sc,A) 



m(A)] , 



(2.10) 

Here, because no small factor of the form 
exp (— ez{A)/x) appears in the prefactor in Ea. (|2.4|) . we 
have assumed that A(A) = Aeg(A) = exp{— ez{A)/x)- 
Then, for simplicity, we have assumed that we can 
bring the latter factor outside of the integral over A in 
Eq. (|2.10| ). and replace ez{A) by a single characteristic 
formation energy ez- Thus we have written 



(A) =exp(-ez/x). 



(2.11) 



We now need to specify the transition rate R{sc,A). 
In earlier work (2ll . |22| | , we wrote this quantity in the 
form 



R(sc,A) = Ro{sc) exp 



(2.12) 



where the function 'D(sc) was equal to the single barrier 
height A when sc = and became vanishingly small at 
large, positive sc- Here, with a range of different A's, 
we must recognize that the stress needed to drive the 
system over a barrier of height A must itself be a growing 
function of A. The simplest such choice is 



1^{sc) = A exp 



ao A 



(2.13) 



where vq, needed for dimensional reasons, can be taken 
to be the same molecular volume introduced in Ea. (|2.7p . 
and flo is a dimensionless number with which we account 
for the uncertainties in the other parameters. If the as- 
sumptions leading to this equation are correct, then ao 
will be roughly equal to unity. The approximation for 
2?(sc) in Eq. (|2.13|) is physically reasonable and mathe- 
matically well behaved except for small A, i.e. for high 
rates, where it will play no role in the analyses to be 
described here 

The prefactor i?o(sc) can, in principle, be any sym- 
metric function of sc- However, since we are considering 



only very small values of sc, this stress dependence is 
irrelevant. On the other hand, because some STZ's are 
complex objects with many internal degrees of freedom, 
we anticipate that the attempt frequency Rq /tq depends 
on both the temperature 9 and the activation barrier A. 
Therefore, in analogy to the STZ creation-rate formula 
in Eq. ()2.8|) . we write Rq = po{9,A). Once again, we 
presume that the activation barriers A are ordinary en- 
ergies roughly comparable to the formation energies ez; 
but that the rates of thermally assisted passage over these 
barriers may be substantially decreased by glassy dynam- 
ics. The A dependence of po is necessary. We expect that 
the transitions over small barriers, near A = 0, are simple 
Arrhenius processes with po — 1- At the other extreme, 
when A is large, we expect po to be small. 

Note that the separation of time scales between the 
configurational and kinetic-vibrational degrees of free- 
dom requires that the internal STZ transitions be no 
slower than configurational rearrangements; that is. 



Pol 



, A) e- >p(0)e- 



(2.14) 



We anticipate that, even if there were a glass transition at 
a nonzero temperature 9o, the prefactor po(^iA) would 
remain nonzero for 9 < 9q although p{9) would vanish. 

These understandings about the rate factors allow us 
to simplify the linearized STZ equations of motion. For 
small Sc, 



C(sc,A)-C(0,A)=po(e,A)e-^/^ 



and 



r(sc,A) -r'(o)sc 



^^0 Sc 
ao9 ' 



(2.15) 



(2.16) 



We stress that C{sc,A) in Ea. (|2.15l) describes only ordi- 
nary thermally activated processes. 
It is convenient to write 



2C(O,A) = 2po(0,A)e 



-A/e 



(2.17) 



and, in most circumstances, to use v as the independent 
variable instead of A. The barrier heights A are energies 
that characterize the configurational subsystem. They 
do not, by themselves, carry dynamic information. In 
contrast, when we write equations in terms of v, we are 
building into them a large amount of information about 
glassy dynamics. Thus, the transformation of variables 
in Eq. (|2.17p is an important feature of the following anal- 
ysis. 
Define 



(2.18) 



Then, up to terms linear in sc, Ea. (|2.10p becomes 



2 To J 



vq Sc 
ao 9 



rh{i') 



(2.19) 
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where m(i/) = m(A). The equation of motion for m, 
Eq. (|2.5p . becomes 



To ifi{iy) = sc - + p) m{v), 



(2.20) 



where we have dropped the term proportional to A on 
the right-hand side of Eq. (|2.5|) . 



III. OSCILLATORY RESPONSE 

A. Prequency-Dependent Modulus G{uj) 

We start by assuming that we are deahng with systems 
that are sufficiently well aged that (A) has approached its 
equilibrium value, and no longer is changing as a function 
of time at rates comparable to experimental oscillation 
periods. We also assume - nontrivially - that the total 
shear rate 7 is simply the sum of elastic and plastic parts: 



(3.1) 



where /i is the shear modulus. 

Denote Fourier transforms as functions of frequency ui 
by m, s, etc.; and let 



Then, 



m(v) = — - 

Similarly, Eq. (|3.1l) becomes 
'iuj Co Wo (A) 



(3.2) 



(3.3) 



where 



p 2 ao To 



[s - jw7)if(w)7] , (3.4) 



J(w) — I dvp{v) V I - 



i w To + p 
Solving for G(iS) — s/^, we find 

AA(w) 



G{lo) 



IUJTq p 



« W To + A J(w) 



where 



and 



A/-(^) = l + !i^^^^(^) + M^AJ(c.) 

P PTO 



2006* 2ao6l ^ 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



To take a first look at these formulas, assume that the 
kinetic viscosity rjK is negligible at low frequencies, and 



use Ea. (l3.6p to compute the Newtonian viscosity associ- 
ated with configurational deformation; 



rjpf = lim 



G(w) 



pro 
A J(0)' 



(3.9) 



It is easy to check that, when the distribution p{iy) is 
sharply peaked at one characteristic value of i/, this for- 
mula reduces to the STZ result in Eq.(5.4) of [2l|. 

Next, return to the expression for G{uj) given in 
Eq. p.6p . and note that this formula cannot naturally be 
expressed as an average over Maxwell modes as in SGR. 
This feature is a result of our kinematic assumption in 
Ea. (l3.ip . plus our assumption that the plastic strain rate 
appearing there is a sum over independent contributions 
from the different kinds of STZ's. 

We can deduce immediately from Ea. (l3.6p . without 
yet knowing much about the distribution function p(i'), 
that the low-frequency structure of G(w) occurs approx- 
imately in the neighborhood of 



Wo 



A J(0) 



To 



_P_ 



(3.10) 



so that Wq is about the same as the viscous relaxation 
rate. However, the structure of the a peak depends sensi- 
tively on the barrier-height distribution, and the approxi- 
mation made in Eq. (|3.10|) cannot replace a full evaluation 
of G{uj). For example, if p(i') were concentrated in a nar- 
row band of values near, say, v = ^ p, there would be 
a relatively narrow peak in G"{uj) at oj = Wq, as given by 
Eq. p.lOp . but G{ll!) would have other structure at and 
above tu — vq/tq. As will be seen, an anomalously broad 
a peak requires that piv) have an appreciable part of its 
support starting well below v — uJaTo- 



B. Barrier-Height Distribution 

At this point, we must begin to discuss the barrier- 
height distribution p(A). Our starting assumption is that 
p(A) is a near-equilibrium feature of the configurational 
subsystem. By this, we mean that it is determined by 
configurational variables, especially the effective temper- 
ature, which themselves may be changing in time while 
the system as a whole moves toward true thermodynamic 
equilibrium. Thus, we start by considering systems that 
are fully aged or are aging slowly, and later will see what 
happens when this assumption fails. 

Since A is measured downward from some reference 
energy, it seems natural to postulate, at least for a range 
of values of A, that p(A) is an equilibrium distribution 
of the form 



p(A) oc e 



+A/A 



(3.11) 



where A is an energy that will be proportional to x 
simple circumstances. For thermosensitive colloidal sus- 
pensions, however, we expect that A will depend on the 
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volume fraction, which itself is a strongly varying func- 
tion of the bath temperature 9. A relation of the form 
A cx X would mean that the distribution of energy bar- 
riers becomes broader as x increases, consistent with the 
fact that X is a measure of disorder as well as energy. In 
the limit of small A and large where po = 1, the factor 
d^jdv in Eq. p.lSp is equal to —O/u, and 



p{v) 



A 



(3.12) 



where A is a normalization factor. 

In the opposite limit of large A and small v, integra- 
bility of p( A) requires that the distribution be cut off for, 
say, A > A*, J/ < i/*. To estimate v* , remember that the 
inequality in Eg. (12. 141) tells us that the STZ transition 
rates must be faster than the rates at which the STZ's 
are created and annihilated by thermal fluctuations. It 
obviously makes no sense to talk about STZ transition 
rates that are slower than the rates at which the STZ's 
themselves are appearing and disappearing. We therefore 
propose that u* be the value of v at which this inequality 
breaks down, i.e. 



2po(e,A*)e- 



2p(6i)e-'=^/^ (3.13) 



Importantly, we do not need to know anything at all 
about po in order to determine v* from the last expression 
on the right-hand side of Eg. p.lSp . This relation pre- 
dicts extremely small values of v* , which are confirmed 
by experiments. The presence of po in the intermediate 
expression simply assures us that we do not have to in- 
voke values of A* that are unphysically larger than ez in 
order to justify small values of v* . 

The cutoff for A > A*, < I/* cannot be infinitely 
sharp; therefore, in this region, we propose to write 

p(A)oce-^/^S (3.14) 
and, equivalently, in the limit of small v, 



(3.15) 



Here, Ai is an as-yet undetermined energy scale, and we 
have assumed (questionably) that the A dependence of 
Po is unimportant in this limit. In general, we expect 
that 1. A much larger Ci - i.e. a much sharper 

cutoff - would require the configurational energy scale 
Ai to be smaller than 6, which seems unphysical. On 
the other hand, substantially smaller values of Ci would 
be inconsistent with our rationale for the choice of v* in 
Eq. p.l3p . (In Sec lIVI we will invoke small values of C,i for 
highly disordered systems far from configurational equi- 
librium, where the preceding rationale for v* is no longer 
valid.) The resulting structure of the function p(A), on 
both sides of its peak at A*, is remarkably similar to 
the distributions deduced from creep measurements by 
Argon and Kuo 
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FIG. 1. (Color online) Experimental data for Vitreloy 4 and 
theoretical comparisons for the storage modulus G'(cj) (red) 
and the loss modulus G"{iu) (blue). The data points are taken 
directly from Fig. 2 of Gauthier et al. Q]- 



Throughout the rest of this paper, we combine the 
high- and low-j/ approximations by writing 



A 



V {{vjv*)'^ + {v* lv)^^\ 



(3.16) 



This is an overly simple, three-parameter representation 
of the barrier-height distribution, with two exponents C 
and Ci determining the large-;/ and small- ^ limits respec- 
tively, and a single crossover value of v* . 

The parameter C controls the high-frequency behavior 
of G{uj). For wTo ^> p, 



G(^) 



1 + A J{u))/iujTo ' 



(3.17) 



and 



J(c.) 



V Piv) dv 



oc 



V dv 



lUJTo Jo V + lOJTo Jj,, U + 

A scaling analysis, using v* <C 2, then tells us that 



(3.18) 



J{io) C(C) 



(3.19) 



where C{Q is a complex constant. Therefore, G"{uj) ~ 
(wTo)~^ in the limit of large uj. 

Suppose, for the moment, that A is simply propor- 
tional to say, A = x/^, where & is a system-dependent 
parameter. Suppose, further, that the system is truly in 
equilibrium, so that x — Then Eg. (13. 121) implies that 
C = b. This parameter "6" , or one nearly equivalent, has 
been computed from first principles in MCT [3l|. Per- 
haps this is a place where we can find a direct relation 
between the STZ and MCT theories. 
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C. Structural and Metallic Glasses 



In comparing these theoretical results with experimen- 
tal data, we look first at the oscillatory response of struc- 
tural and metallic glasses, for which tq is of the order of 
picoseconds, lutq ^ 1, and rjK is negligible. The inter- 
esting behavior occurs at temperatures near or slightly 
above the glass temperature. 

Our principal sources of information about the oscilla- 
tory responses of structural and metallic glasses are the 
papers by Gauthier et al., in particular [2|. These au- 
thors show that the functions G{uj), for a wide variety 
of noncrystalline materials at their glass temperatures, 
have very similar behaviors. Specifically, the loss modu- 
lus G"{lu) has a broad peak at cJq, and drops off at high 
frequencies like uj~'^ as predicted in Eg. (13.191) . For metal- 
lic glasses, Gauthier et al. find ( = 0.4. 

In FigUl we show G'(a;)//x and G"{oj)/ii aspredicted 
by Eg. p. 61) . along with data from Fig. 2 of for the 
metallic glass Vitreloy 4 at its glass temperature Tg. 
In estimating the theoretical parameters, we have used 
Tg ^ 600 i^, Tq 2 X 10^12 ggc.^ and ^ 50GPa. 
To make approximations for the other parameters in 
Eg. (l3.6p . we note that, if the volume vq is of the or- 
der of a few cubic nanometers, then the ratio vofi/Og, 
and thus the prefactor in Eq. (|3.8|) . is approximately 10^. 
Then, to estimate {A{9g)) ~ exp{— ez/dg), we assume 
that 9g is the same as the steady-state value of the ef- 
fective temperature, usually denoted by xo, for systems 
driven persistently at shear rates much slower than Tq^. 
The ratio Xo/^z may be a universal quantity; it usually 
turns out to be in the range 0.1 — 0.2. We therefore es- 
timate Og/ez ^ 0.15, implying that {A{9g)) ~ 10~^, and 
therefore that A ^ 10. As a result, Eq. p.l3|) tells us that 
ly* = 10-3 

The theoretical curves in Fig[l] have been computed 
using Ci = 1, pi9g)/To = 1.25 x lO'^ sec.-\ and A = 25. 
In effect, we have set eo/oo ~ 2.5 in Eg. p.Sp . which is 
well within our theoretical uncertainty. Note that these 
parameters imply that v* ~ 10"^^, which, as predicted, 
is extremely small compared, for example, to its upper 
limit ai v — 2. So far as we can tell from numerical explo- 
ration, this small value of v* is sharply determined by the 
experimental data. Because p{u) is varying very rapidly 
near v — v* , changing v* by even a factor of 2 ruins 
the fit to the data; and the curves become qualitatively 
wrong if v* is changed by an order of magnitude in ei- 
ther direction. Our theoretical curve for G"{uj) in Fig.l 
is somewhat less smooth than the experimental curve, 
which may be a result of our overly simple form for p^v) 
in Eg. (l3.16p . The other discrepancy is that the experi- 
mental storage modulus G'{uj) is approximately linear in 
w at low frequencies, instead of being proportional to 
as predicted by our theory as well as by Maxwell models. 



100 




0.01 0.1 1 10 100 1000 10" 

cj [rad/sec] 




o) [rad/sec] 

FIG. 2. (Color online) Experimental data and theoretical 
comparisons for the storage modulus G' (uj) (red) and the loss 
modulus G"{uj) (blue), for three different suspensions of ther- 
mosensitive particles, as reported in 9]. The values of the 
parameters are listed in the text. 



D. Colloidal Suspensions 

The rhcology of soft, colloidal suspensions differs from 
that of structural and metallic glasses in at least two 
important respects. First, in colloidal systems, the ap- 
proach to jamming near a glass transition is controlled 
more sensitively by the volume fraction than by the tem- 
perature. Second, the microscopically short molecular 
vibration period in structural glasses is replaced in col- 
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loids by the very much longer time scale for Brownian 
motion of the particles. As a result, the high-frequency 
cutoff at = 2 is probed in rheological experiments, and 
the kinetic viscosity fixi'^) is relevant at accessibly high 
values of w. When samples are prepared by subjecting 
them to strong shear stresses, their shear rates become 
comparable to . Therefore, they become highly dis- 
ordered and their effective temperatures x can become 
arbitrarily large (32j . 

In order to use Ea. p.6p to evaluate G{lS) for colloidal 
suspensions, we need an expression for the frequency- 
dependent kinetic viscosity f}K{uj). Here, we follow Lion- 
berger and Russel [s^, who show that G'(a;) « w^/^ in 
the limit of very large w. Their analysis is based on the 
idea that, in order to satisfy the no-flow boundary condi- 
tion at the surface of a colloidal particle, the surrounding 
fluid must form a diffusive boundary layer whose thick- 
ness scales as The viscous force on this particle 
therefore decreases like ryx (w) ^ lj"^/^, and the resulting 



stress is i w r?_ff (i^) 



1/2 



We regularize this expression 



at low frequencies by assuming that 



(3.20) 



where c is a dimensionless constant, and tk is a viscous 
time scale. 

In Figl^J we show three examples of how the STZ the- 
ory developed here is capable of reproducing the experi- 
mental results of Siebenburger et al. [9]. These authors 
explored a range of effective volume fractions 0e// and a 
wide range of frequencies a; (as well as steady shear rates 
not discussed here) by using suspensions of thermosensi- 
tive particles (polystyrene cores with attached networks 
of thermosensitive isopropylacrylamide molecules). They 
have used the effective size of the particles and the sol- 
vent viscosity to estimate that the Brownian time scale 
is To - .003 sees. [H]. 

The sequence of increasing volume fractions, from 
the top panel to the bottom in Figl^l is 0e// = 
0.518, 0.600, and 0.626. The theoretical parameters, de- 
duced by fitting the data and listed in the same order, 
are p = 0.04, 3 x lO"'^, lO'^; v* = 0.001, lO'^p, 10"^; 
C = 1.0, 0.5, 0.38; A = 200, 40, 17; = 12, 35, 45 Pa; 
and TK = 0.004, 0.002, 0.002 sec. In all cases, Ci = 1 
and c = 0.1. 

The trends are interesting. The sequence of examples 
starts in the top panel of Figl2] with a system whose 
relatively small volume fraction puts it well away from 
the glass transition. It is effectively a liquid; p = 0.04 
means that there is relatively little super- Arrhenius sup- 
pression of the structural relaxation rate. As the systems 
become more glassy in the middle and bottom panels, p 
decreases rapidly. C also decreases, as if the temperature 
9 in Eg. p. 121) were decreasing; but it is the increasing 
volume fraction that must be causing the energy scale 
A to increase. At the same time, the shear modulus /i 
increases and the viscous time scale decreases slightly 
as the systems become stiffer. A is very large for the liq- 



uidlike example in the top panel, implying that the STZ 
density is large in this system. On the other hand, the 
values of A, and the relation between v* and p for the two 
nearly glassy cases, are almost the same as the analogous 
estimates for bulk metallic glass in the preceding subsec- 
tion - despite the fact that the underlying time scales for 
these systems differ by nine orders of magnitude. 



IV. STRAIN RECOVERY AND AGING 
A. Strain Recovery: Basic Aspects 

The success of the STZ theory in accounting for a wide 
range of oscillatory measurements gives us some confi- 
dence that our fundamental concepts are correct. We 
turn now to a class of experiments whose interpretation 
is not nearly so simple. 

Although strain-recovery experiments technically 
probe the linear response theory outlined in the preceding 
sections, they are qualitatively different from the oscilla- 
tory experiments because the mode of sample prepara- 
tion is necessarily nonlinear. In order to induce a strain 
whose recovery can be observed, these systems are de- 
formed at stresses well above the yield stress. As a result, 
we cannot use near-equilibrium concepts to estimate the 
initial barrier-height distribution as we did for the oscil- 
latory analysis. Nor can we assume that the initial value 
of the effective temperature x is close enough to a steady- 
state value that we can neglect its time dependence. 

Strain recover y ni easurements have been carried out by 
Belyavsky et al. Il3l| for metallic glasses, and by Purnomo 
et al. ji] for colloidal suspensions. We focus on the more 
recent colloidal experiments described in Q ; and we fur- 
ther specialize to results for the thermosensitive poly-N- 
isopropylacrylamide microgel suspension that Purnomo 
et al. call "P-1". A representative data set taken from 
Fig. 4b in Q is shown here in FigEl 

In both metallic glass and colloidal cases, the systems 
first were deformed at high stress, and then unloaded 
quickly. After an initial, almost instantaneous, elastic 
relaxation, the plastic strain decreased slowly at zero ap- 
plied stress, as seen here in Fig|31 We compute this plas- 
tic relaxation by setting sc = in Eq. (|2.20p . obtaining 

To m{h') ~ —{v + p) m{h'); fh(y) = rho e~ ('^+'') */'^o ^ 

(4.1) 

where, for simplicity, we assume that the initial values of 
all the 771(1/) are equal to ttiq. Then, using Eqs. (j2.19[) and 
p.ip . in both cases with sc — and no kinetic viscosity, 
we have 

to , . > / 
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-—2- (A) / dvp{v)vm{v,t) 



= _£^(A)/ diyp{i/)i/e~^''+P^'/^\ (4.2) 
2 To Jo 

A first, rough comparison between these formulas and 
the experimental strain-recovery data tells us that we are 
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FIG. 3. (Color online) Strain recovery as a function of time. 
The blue data points are taken directly from Fig. 4b of [^. 
The strain is measured from its value after unloading and 
after an initial, almost instantaneous, elastic relaxation. The 
red theoretical curve has been computed using the parameter 
values discussed in the text. 



seeing an entirely different distribution p(i^) than the one 
that we used to interpret the oscillatory data in Sec lIIII 

One clue about what might be ha ppe ning is contained 
in a paper by Rodney and Schuh [35|- These authors 
used Monte Carlo simulations to track the distributions 
of thermally activated events in a binary, two dimen- 
sional, glassy material. When their system was allowed to 
age at sub-yield stresses, they found decreasing numbers 
of increasingly high activation barriers. In contrast, when 
they persistently deformed their system at high stresses, 
they found larger numbers of lower barriers. Karmakar et 
al. [3^ have reported similar observations about barriers 
seen in athermal quasi-static simulations. 

In the language used here in Scc lIIIi with the notation 
introduced in Eq. p.l6|) for the distribution p{i') , the first 
case discussed by Rodney and Schuh resembles the 
small-i^* distribution at small x that we used to describe 
the oscillatory responses of well aged systems. In their 
second case, they drove their system to large x, i.e. to 
high energy and high disorder, where they saw large num- 
bers of low activation barriers. This case corresponds to 
a large value of v*, with the weight of the distribution 
in the range ly < ly* , and with p{iy) ^ The ap- 

proximation for v* in Eq. (j3.13p is not relevant in this 
far-from-equilibrium situation, where the barrier heights 
near v* are too small to be limited by thermally induced 
configurational fluctuations. 

Suppose, for the moment, that we can ignore the time 
dependence of (A) in Ea. (l4.2p . Then 
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FIG. 4. (Color online) Strain increments induced by small 
stress steps at waiting times of = 30, 600, 3, 000 and 10, 000 
sees, indicated respectively, from top to bottom on the left- 
hand side, by blue, green, red, and black data points and 
theoretical curves . The data points are taken directly from 
Fig. 5a of ,8]. The parameters used in computing the theoret- 
ical curves are discussed in the text. 



Thus, after unloading to zero stress, the strain ultimately 
relaxes by an amount 700, and does so at a rate that is 
no slower than exp(— pt/ro). For large values of ^/tq, 
the integral over v in Eq. (j4.3p is dominated by the part 
oi p{v) for v < v* , and therefore: 



lit) « -7oo + ^ (A)i,e-''*/-° / 

Jo 



-7o 



2 

const. 



ro/t 



- pt/To 



[t/ro] 



(4.4) 



(4.3) 



If Ci is small, that is, if the distribution p(z/) is nearly 
flat down to small values of then strain recovery will 
appear to be logarithmic out to times of the order of tq /p, 
which is the behavior described by Belyavsky et al.|l3|. 



B. Stress-Step Experiments and Structural Aging 

The strain-recovery data for the colloidal suspensions 
of Purnomo et al. [8] can be fit reasonably well by us- 
ing Ea. (l4.4p with values of in the range 0.2 — 0.3. 
With only this information, we might have concluded 
that these authors were seeing just the rh relaxation de- 
scribed by Eq. (j4.1l) . and not true structural aging in 
which STZ's are created and annihilated. But these au- 
thors went further in important ways. After preparing 
their samples by rapid straining and abrupt unloading, 
they allowed the strain to relax at zero stress for various 
waiting times tw, and then measured the response to a 
small step in the stress whose magnitude dsc was much 
less than the yield stress. Their results are shown here 
in FiglH 
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At least two aspects of these results can be understood 
only in terms of true structural aging. First, the initial, 
instantaneous, elastic increase in the strain, following the 
stress step, decreased with increasing aging time, indicat- 
ing that the shear modulus fi was larger in the older sam- 
ples. The shear modulus must be a decreasing function of 
X (just as it is a decreasing function of the ordinary tem- 
perature 9). Thus, we conclude that x ~ an intrinsically 
structural variable - decreased during aging. Second, 
Purnomo et al. found that, during the initial responses 
to the stress steps, the functions j(t) for different waiting 
times could be scaled approximately onto each other if 
they were plotted as functions of t/t^, where i = is the 
time at which the stress step was applied. Both of these 
features can be seen here (roughly) in FigU) 

To interpet these behaviors, consider the equation of 
motion for 7 for times t > 0. If Eq. (|4.ip were accurate 
throughout the aging process - which we will argue is not 
exactly the case - then 



m{i', t — 0) ~ rho e 



-(1'+P)t™/T0 



(4.5) 



We insert this formula into Ea. (|2.19p . and find that the 
subsequent rate of deformation is 
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{A{t + t^)) 



To 
P 



di'p{v) V |- 



TO e 



{v+pYItq 



V + p V + p 



where 



'>n = —niQ, 



£0 vq Ssc 
2aae 



(4.6) 



(4.7) 



The first term inside the brackets, proportional to fh, is 
the continuation of the original strain-recovery formula 
following the waiting time i^. Of the two terms pro- 
portional to the dimensionless stress step cr, the first is 
the viscous response, consistent with Eq. (j3.9p . which is 
nonzero so long as p is nonzero. The last term tells us 
that, when p is small, this system initially behaves like 
a conventional yield-stress material; that is, for small 
stresses, it fairly quickly becomes jammed with a step- 
wise, fixed increase in the plastic deformation. In con- 
trast, the viscous strain proportional to p continues to 
grow linearly in time, and thus can become arbitrarily 
large. 

A crucial point is that each of the terms in Ea. (l4.6p 
is proportional to the aging factor {A{t + t^)) — 
exp [— ez/xi^ '^t'w)]- comparisons with experimental 
data, we must compute this function exactly by solv- 
ing Eq. (|2.9p numerically; but it is useful to look first at 
an analytic approximation. For simplicity, assume that 
CA = ez , and also assume that we are interested primar- 
ily in situations where x, while decaying, is still much 
larger than 9. Then, in the limit if < » 1, Ea.(|2^ has 
the asymptotic solution 



p-ez/x(t) 



1 



Kt \xi{Kty 



K 



nez P 
To 9 ■ 



(4.8) 



Therefore, apart from the logarithmic correction in 
Eq. ()4.8p . all the terms on the right-hand side of Eq. ()4.6p 
contain the factor [t + t^)^^ k t^^ for t t^, which is 
roughly consistent with the observed behavior. 



C. Preliminary Comparisons between Theory and 
Experiment 

We turn, finally, to quantitative comparisons between 
our theory and the experimental data. We will argue 
that, although we are able to achieve fairly good agree- 
ment with the experiments, the remaining discrepancies 
indicate that the theory is missing an important ingredi- 
ent. Specifically, we still need an equation of motion to 
describe the time dependence of the barrier-height distri- 
bution p(v) during structural aging. 

We have used Eq.(I12]), along with Eqs.(I13) and (gH) 
for strain recovery and stress-step response respectively, 
to interpret the experimental data shown in Figs. [3] and 
m Figure [3] shows data for strain recovery taken directly 
from Fig. 4b of l8|. Figure 2] shows data taken directly 
from Fig. 5a of for responses to stress steps applied 
after waiting times of 30, 600, 3, 000, and 10, 000 sec. 

Consider first the comparison between theory and ex- 
periment in Figl21 Using values for the particle radius 
and solvent viscosity given in 's'l , we have estimated that 
To = 0.1 sec. In solving Eq. (|2.9p for the ratio x(i)/ez, 
we have set ca — ez and, assuming that the system is 
near its glass temperature, have used 9/ez = 0.15 as dis- 
cussed in Sec. IIIII C. We have evaluated the rate factor K 
defined in Ea. (j4.8p by assuming that the shear modulus 
p{x) is a linear function of the form p,o — pi x, and then 
fitting this relation, with x — xi^w), to the values of the 
initial elastic increments seen in Fig. 5a of @. The re- 
sult is that K ~ 10~^ sec.~^, which already is interesting 
because it implies that the rate at which x relaxes to 9 
in Ea. (j2.9p is comparable to the rate of strain recovery 
seen in Fig. [31 With this value of K, and with the pre- 
ceding estimates of tq and 9/ez, we guess that k 0.1 
and thus have chosen p ~ 10~^. We also - by necessity 
- have assumed that the volume fraction remained fixed 
throughout these experiments. 

The theoretical curve shown in FiglS] has been com- 
puted using p = 3 X 10"^, K = 10^^, fh = eofaojl — 
0.157, C = 0.4, Ci = 0.32, x(0)/ez = 1-0, and v* = 0.5. 
With these choices of parameters, the agreement between 
theory and experiment seems good; but, in fact, it is quite 
problematic. Most obviously, i/* = 0.5 is vastly different 
from the value v* ~ 10^'^p ^ 10"^ that we expect to 
find in nearly equilibrated systems. A related problem 
is that we have assumed that the barrier-height distribu- 
tion p{y) remains fixed during this process; but p{y) is a 
structural property, and must change during structural 
aging of the kind implied by decreasing values of x(^)- 
Aging somehow causes the system to make a transfor- 
mation from its initial, highly disordered state, with a 
barrier-height distribution characterized by large values 



11 



of X Etnd V* , to more ordered states with substantially 
smaller values of those parameters, especially the latter. 
As yet, we have no theory of how this strongly nonequi- 
librium process occurs. We do not know, for example, 
whether v* might be a near-equilibrium function of the 
instantaneous value of x, or whether it relaxes toward its 
near-equilibrium value on some dynamically interesting 
time scale as happens in [2Q| , or whether any such weakly 
nonequilibrium mechanism is operative. 

Turn now to the stress-step data in FigH) For well aged 
systems such as those for = 3, 000 and 10, 000 sees., it 
seems safe to assume that the configurational subsystems 
have settled into states that are much the same as those 
discussed in Sees. |llI]C and D, with v* ~ lQ~^p. Because 
of this qualitative change in p(v), there is no reason to 
expect that the system's orientational memory, carried 
by the parameter m in Ea. (j4.6|) . retains its original value; 
thus we use it as a fitting parameter, and find that it does 
undergo modest changes during aging. 

We have used the following procedure to determine the 
theoretical parameters for the four cases shown in FigU) 
As we did for the strain-recovery curve in Fig|3l we have 
computed each of these curves separately by choosing 
constant values of p, K, i/*, C, Ci, and m, rather than 
trying to guess how these parameters might depend on 
the time. We have kept p — 3 x 10^^ and K ~ 10^^ 
from the strain-recovery analysis. Despite the fact that 
v* changes dramatically during the early stages of aging, 
we have chosen i/* = 10~^ p for all four cases. This re- 
lation has a theoretical rationale, and the data does not 
compel us to choose otherwise - not even for = 30 
sec. We looked first at the most fully aged system, with 
= 10,000 sec. Here, as argued in Sec lIIII D. we have 
assumed that (i = 1. We then have fit the data by choos- 
ing C = 0.21, m = 0.09, and a = 0.07. This value of a 
must be common to all four curves in FigU) 

For the three other waiting times, the values of the pa- 
rameters are: for t^i, = 3, 000 sec, C = 0.29, Ci = 1, "i, = 
0.092; for ty, = 600 sec, ( = 0.28, Ci = 0.5, fh = 0.125; 
for = 30 sec, C = 0.23, Ci = 0.2, m = 0.130. As 
expected, the values of m decrease with increasing t^^ 
but they differ from the value (0.157) deduced from the 
strain-recovery data by less than 50%. The numerical 
results for = 30 and 600 sec. are sensitive to the value 
of Ci , which seems to increase with increasing i^, . We see 
no pattern in the values of which may be artifacts of 
our parameterization of ^{v) or of our parameter-fitting 
procedure. They might also be explained by changes in 
the volume fraction during aging. There are other exper- 
imental uncertainties. The data shown in Fig. 5a of [8| is 
very noisy. As a result, our estimates of the initial elastic 
strains in Fig|4] are uncertain; and we do not know how 
seriously to take other discrepancies between theory and 
experiment at later times. 

One systematic theoretical error is that, for the two 
cases with shorter aging times, the rising viscous strain 
is too large at the latest times observed. We have em- 
phasized this behavior in Fig|3] by extending the time 



axis an extra decade beyond the data. The formula for 
the viscosity in Ea. (l3.9p depends on ^{v). In these cases, 
p{v) has been chosen primarily to fit the earlier stages 
of the measurements; thus, it is not surprising that the 
viscosity is wrong at the later stages. 

The most interesting feature of these results is the huge 
change in v* that apparently occurs in the first ten sec- 
onds or so after unloading. It appears that the popula- 
tion of STZ's with low activation barriers in the strongly 
disordered inital state disappears quickly - much more 
quickly than can be described by near-equilibrium x dy- 
namics - so that these STZ's are not available to respond 
to the only slightly delayed step in the stress at = 30 
sec. This observation suggests that, even if we could 
deduce an equation of motion for v* , we might not cap- 
ture essential features of the configurational dynamics in 
this regime. We then must ask, if p{v) is such a strong 
function of aging time, why are we able to describe strain 
recovery with a time- independent p(y)l Perhaps the be- 
havior seen at long times in Figl3]is sensitive only to the 
population of slow STZ's, and perhaps this population 
remains unchanged by aging during the times of interest. 
But such an explanation does not yet provide a basis for 
developing a predictive theory of these phenomena under 
far-from-equilibrium conditions. 



V. SUMMARY AND CONCLUSIONS 

The analyses described here lead us to a number of 
conclusions and assertions. Some of these, such as (3) 
and (5), are shared by SGR and MCT, but emerge here 
in different ways. Others seem to differ more strongly 
from earlier formulations. 

(1) Linear response in glassy systems is determined by 
thermally activated processes. Long time scales can be 
explained, in analogy to structural relaxation in glasses, 
by super- Arrhenius reduction of attempt frequencies as 
opposed to anomalously high activation barriers. 

(2) STZ dynamics and kinematics, with a nontrivial 
barrier-height distribution, produces a non-Maxwellian 
form of the frequency-dependent modulus G(lS). 

(3) The loss modulus G" (iS) automatically has an a 
peak in the neighborhood of the viscous relaxation 
rate, independent of any specific distribution of barrier 
heights. 

(4) For fully aged or slowly aging systems, the principal 
features of the STZ barrier-height distribution j)(y) can 
be determined by assuming that the system is nearly in 
equilibrium with the current values of various configu- 
rational variables, most importantly, with the effective 
temperature x- However, a detailed, first-principles 
relation between piv) and x remains undetermined 
by the analysis presented here, even for equilibrium 
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situations. 

(5) Our near-equilibrium analysis of the barrier-height 
distribution in well aged glassy materials implies that 
any such material exhibits an anomalously broad a 
peak near its glass temperature. Metallic and structural 
glasses, and colloidal suspensions, all share remarkably 
similar properties in this regard, despite the enormous 
differences in their underlying time scales and internal 
dynamics. 

(6) Some, but not all, of the dynamics of slow processes 
such as strain recovery can be attributed to internal 
orientational relaxation, rather than to true structural 
aging. 

(7) True structural aging is characterized by rear- 
rangements of configurational degrees of freedom. In 
particular, the barrier-height distribution must change 
during structural aging or persistent deformation. At 
present, we still need an equation of motion for p{i'), 
perhaps analogous to the SGR equation of motion for 
the distribution over barrier heights and local strains. 
The nonequilibrium dynamics of the barrier-height 
distribution is the most important of the issues raised 
and left unresolved in this paper. We hope to return to 
it in the near future. 

The analyses presented here have not touched on sev- 
eral other topics that often are included in rheological 
discussions. For example, we have not discussed stress 
relaxation in fixed-strain experiments. Among the sub- 
tle issues that arise in this connection is that the Fourier 
transform of the stress relaxation function is not neces- 
sarily the same as the function G{lo) defined here. The 
difference is that our systems contain relevant internal 



variables such as rh{y) and x, which relax on the same 
time scales as the stress, so that conventional theories 
with history-independent relaxation kernels are inappli- 
cable. 

We also have not discussed nonlinear plasticity, be- 
cause that topic has been the main focus of recent STZ 
theories. For example, the curves of stress versus strain 
rate shown by Siebenburger et al 'gj] for colloidal suspen- 
sions are qualitatively similar to those that we have dis- 
cussed in earlier papers on bulk metallic glasses [2l|, [l^ . 
We believe that we understand those phenomena, and 
that they are not especially sensitive to the barrier-height 
distributions. However, we will need to invoke nonlinear 
mechanisms to understand the sample-preparation meth- 
ods used by Purnomo et al fs'] . 

Finally, we have considered only systems in which the 
ambient thermal noise is no n- negligible. We have not yet 
extended these ideas to strictly athermal systems such as 
foams or granular materials. It will be interesting to do 
so. 
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